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We present a semi-analytical approach to the interaction of two (originally) Kerr black holes 
through a head-on collision process. An expression for the rate of emission of gravitational radiation 
■ is derived from an exact solution to the Einstein's field equations. The total amount of gravitational 

radiation emitted in the process is calculated and compared to current numerical investigations. We 
find that the spin-spin interaction increases the emission of gravitational wave energy up to 0.2% of 
the total rest mass. We discuss also the possibility of spin-exchange between the holes. 
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The two body problem is one the most important and difficult challenges in general relativity. As a particular 
case, the head-on collision of two black holes has recently gained new insights from several techniques, namely, 
numerical 0] , perturbative || , post-Newtonian H and semi-analytical methods [Q . These in turn concentrated their 
initial efforts in studying the collision of two Schwarzschild black holes. A semi-analytical approach developed Q to 
t-H ■ treat the Schwarzschild case was based on an exact solution of Einstein's field equations. A particular case of this 
solution describes a geometry that can be interpreted as a static axisymmetric spacetime with two black holes plus a 
conical singularity between them. It is then possible to obtain the "force" of attraction between the holes and their 
"acceleration". These expressions were used to calculate the rate of emission and the total amount of gravitational 
wave energy released in the process. The remarkable agreement between these results and those of the numerical 
treatment seems to indicate that the static acceleration contribution is very important for the head-on collision 
!h ■ problem. This conjecture is emphasized by the result that, in the context of perturbation theory, the gravitational 

wave emission of boosted perturbed black holes is very similar to that of the static case B . 
^ , The previous success with the semi-analytical approach strongly motivated the present development of a similar 
treatment of a toy model for the head-on collision problem of two Kerr black holes. We start with an exact solution 
of Einstein's field equations that can be interpreted as a stationary axisymmetric spacetime with two collinear Kerr 
black holes and a singularity on the axis of symmetry. This solution has some free parameters that determine the 
nature of the singularity. This singularity has been shown to be either a spinning rod between the holes or a pair of 
spinning strings extending from each hole to infinity, or both in the most general case ||. For each of these cases, one 
can obtain a stationary "force" of attraction and the corresponding "acceleration" from which we calculate the rate of 
emission and the total gravitational wave energy released in the collision when we "cut" the strut or the strings and 
let the holes collide. It is important to point out that in the Kerr case there is also a spin-spin interaction between 
the holes || , which is always a repulsive contribution to the "attractive force" . 

The plan of this paper is as follows: In section II we discuss the "stationary" force between the holes with axisym- 
metric spins. In section III the equations of motion for the holes are obtained. In section IV we derive an expression 
for the rate of emission of gravitational radiation and calculate the total output of gravitational wave energy released 
for several initial configurations. These results are summarized in one table. Section V is dedicated to the possibility 
of the holes' spin variation during the collision process and section VI to our conclusions. 

II. THE TWO BODY INTERACTION FORCE 

The line element of a stationary axisymmetric spacetime is given by the Lewis metric J^] 

ds 2 = exp(2V>) (dt - w defy 2 - exp(-2^) [exp(» {dr 2 + dz 2 ) + r 2 dq9] (1) 
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in cilindrical coordinates < r < oo, — oo < z < oo and < 4> < 2ir . The metric functions ip = il>{r, z), v — v{r, z) 
and w = w(r, z) satisfy the vacuum Einstein's field equations: 



Ar + A/r + tf„ - {wl + w 2 z ) , (2) 

w rr -w r /r + w zz = -j(w r ^ r + w^ z ) , (3) 

,V=r(#-tf) , (4) 
^ = 2np r ipg - ° xP 2r ~ ^ ■ ( 5 ) 

The metric function ^ is obtained, up to a constant, by quadrature once the solutions tp and w for the non-linear 
coupled system of equations (|^) and (^) are known. These equations are also the integrability conditions for the 



rdrd<pdz [M 



equations (^) and (^) and they can also be obtained from the action J + ipl) — cx ^^ 

Asymptotically ip plays the role of a Newtonian potential and w/A of a body's angular momentum per radial distance. 
w is related to the usual twist potential x by r \ r = w z and r x z — —w r . 

Solutions of the system of equations @-(|^) have been obtained using Backlund transformation and inverse scat- 
tering techniques In general these solutions have singularities. We point out that the choice of the integration 
constants in the solution is a crucial step towards the interpretation of the nature of the singularity and consequently 
of the resulting spacetime geometry. In this paper we are going to deal with two possible configurations. One that 
corresponds to two collinear Kerr black holes with a twisted spinning strut between them and the other case, that of 
a couple of spinning strings extending from each of the holes to infinity. 

The strut and string's energy momentum vanishes everywhere except at the axis. Furthermore, their effective mass 
is null and their angular momentum per unit length is given by wexp( , 0)/4 at r = p"2] ] Q. Let mj , ra 2 be the 
masses and a\, a 2 be the angular momenta per unit mass of the holes (hereafter called simply spins), respectively. It 
has shown that the compression force on the strut is given by Q 

p m im2 , . 

" strut — , 2 ( i \2 i I \1 ' \ ) 

d z — (mi + m2) + (ai — 0,2) 

where d is the coordinate distance between the center of the holes. We assume that this is also the "stationary force 
of attraction" between the holes if the strut is removed. In order to get rid off the semi-infinity strings the spacetime 
has the constraint ^ = ^ so it is applicable to the parallel spin case only. 

The anti-parallel spin case can be treated using the tension force on the string as the "stationary force of attraction" 
between the holes if both strings are cut. The force of attraction is given by 

iWngs- ^ _ ^ _ m2)2 + (oi _ a2)2 (<) 

which is opposit to the tension on the strings. 

Both forces are singular at a finite separation d , nevertheless our computation stops before that. Recall that 
(Tfc = \J ml — a? , TOfc > afc , k = 1, 2 is the coordinate distance between the center of each isolated hole and its infinite 
redshift surface along the axis. In this paper we consider d 2 > (mi +m 2 ) 2 — (ai — a 2 ) 2 > (<7i + (T2) 2 . Recall also that 
these coordinates distances correspond asymptotically to the distances measured by an observer at infinity. 

Note that the spin-spin interaction reduces the attraction intensity unless the holes are co- rotating (a\ =02). If the 
holes have anti-parallel spins ( a\ = —02), the spin-spin repulsion has a maximum |g]. But it gives always a repulsive 
contribution to the force. For large distances the spin-spin repulsion is m\m2{ai — a 2 ) 2 /d 4 ■ 



III. EQUATIONS OF MOTION 

We now investigate the time evolution of a two-body system interacting through the force (^|) applicable only to 
the parallel spins black holes. We assume that Newton's equations of motion are adequate to describe the black holes 
motion and that their spins and the masses are constants. We then have a Newtonian two-body problem that can be 
transformed into a central force problem with a fixed center of mass. Let M = mi + 1112 and /1 = mim^/M be the 
total and reduced mass of the system, respectively. Let = \<Xi — 02 1 be the relative spin of the holes. Their relative 
acceleration is given by 
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z 2 -M 2 



(8) 



Let z — z/M , a = clr/M and t 
a < 1. 

The first integral of equation (p 



t/M. Recall that the spin of a hole is constrained to be less than its mass, so 



\z +V 



(9) 



where e is the total energy per unit of reduced mass, z = dz/dt and the potential V is 



V{z,a) = ~]n 



(10) 



where a = Vl — a? . 

We now investigate the time evolution of a two-body system interacting through the force given in ([?]) applicable 
to the anti-parallel spins black holes. We remark that this space-time is not asymptotically flat because of the semi- 
infinity strings. Even so we assume that Newton's equations of motion are adequate to describe the black holes motion 
and that the spins and the masses of the holes are constants. 

Their relative acceleration is given by 



a z 



mi + f«2 



z 2 — (mi — m 2 ) + a R 



Let z = z/M , a — cir/M , 5 = (mi — 7712) /M and t — * t/M. The first integral of equation ( |ll| ) is similar to 
the potential 



(11) 
with 



V(z, a, 6) 



z-/3 

arctan 
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for S 2 > a 2 
for 5 2 = a 2 
for 5 2 < a 2 



(12) 



where (3 = V ' 5 2 — a 2 . 

Therefore we can treat both, parallel and antiparallel spins cases, using the potential (|l^) since a = f3 if we set 
5 = 1 . If S 2 = a 2 the potential becomes the Newtonian. 



IV. EMISSION OF GRAVITATIONAL RADIATION 



For the computation of the gravitational radiation luminosity (rate of emission of gravitational wave energy) we use 
the standard quadrupole formula . In our model this is a good first approximation to get an estimate of the total 
amount of gravitational radiation during the head-on collision. The formula takes an average over several wave-lenghts 
and we integrate the luminosity along the whole collision path up to a point where the holes coalesce. 

For a two particle system separated by a distance "z the luminosity is given by 



dEg 
dt 
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(13) 



Our main concern is the energy emitted along the collision process. Therefore we need the rate of emission of 
gravitational wave energy per unit length. For a head-on collision of two bodies with equation of motion (0) the 
quadrupole formula (|l^) yields 



dE g _ dE g 1 



dz 



dt 
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(z 2 -P 2 ) 

2 (z 2 +m\ 2 ) 
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for S 2 > a 2 , 
for S 2 = a 2 , 
for 5 2 < a 2 . 



(14) 
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Note that, for S 2 > a 2 , as z is decreasing in the collision process,-^ reaches a maximum at z w , then it 
decreases to zero at z = V3/3 and finally it diverges at z = /3 . The behavior of the rate of emission of gravitational 
radiation as a function of the distance between the holes is similar to the Schwarzschild case PJ . The main difference 
is the fact that the two Kerr black holes horizons touch each other at ~z = (<j\ + er 2 ) /M < (3 < a and at this point 



dE 



—^M- is finite, whereas it diverges at the horizon's touching point in the Schwarzschild case. In any case, as before, 
the description of the head-on collision problem by this approach is valid only for z > Zfinal = 0/ tanh (§j> the point 
where the holes would reach the speed of light if released at the rest from infinity. 

For the S 2 < a 2 case, as ~z is decreasing reaches a maximum at ~z — and it is finite everywhere. Again the 

description is valid only for z > Zfi na i = \P\ / tan (^Pj due to the speed limit. 
The total amount of gravitational wave energy emitted is then given by 



AEg (Zjnjtial, Zflnal) 



dz 



-dz 



where ^initial is the point of release from rest. We make a restriction to the cases — 1/2 < e < , which corresponds 
to Zfinai < ^initial < oo. Our results are summarized in the table below, where the total amount of gravitational wave 
energy is given as a function of e and (3 in units of x /i 2 /Af. 

The efficiency of the emission of gravitational waves decreases with j3 for S 2 > a 2 , increases with 
and decreases with e . Note that the anti-parallel case can emit as much as 3% of the rest energy fi 2 /AI . We see for 
example that the collision of equal masses extremal (m = a) black holes with anti-parallel spins (f3 — i) can release 
at least 3 times more gravitational wave energy than the similar collision with parallel spins (/? = 1). 



for S 2 < a 2 



V. TIME VARIATION OF SPINS 



Let us now turn to the possibility of time variation of the holes' spins. Since both strut or strings attached to the 
holes have angular momentum per unit length, they may affect the angular momentum of the holes along the collision 
path. It would be a mutual Lense-Thiring effect or the dragging of each other's inertial frame. The gravitational 
radiation does not carry out angular momentum in the quadrupole approximation due to the axisymetry jl . Therefore 
it is reasonable to expect angular momentum conservation in the configuration. 

In the parallel case the strut's total angular momentum is given by B 



Jsi 



2m 2 di ( d + m 2 + mi) (d — a% — ct 2 ) 
d 2 — {mi — m 2 ) 2 + (ai — a 2 ) 2 



has been used. Let the total angular momentum be I = J s) 



-raidi +m 2 a 2 . Let us assume 



The constraint — = 

that I is conserved during the collision process. Then a% and a 2 have to change with the separation distance d. Let 
a = 01/mi = a 2 /m 2 , S = m 2 — mi and set mi + m 2 to one. In order to keep £ constant, a must be a real root of 



(5 2 + 2zS 2 



1 



(Z+l) (i-s 2 ) 



8 2 (l + 6 2 ) 
(z + l){l-S 2 ) 



(1 - a 2 ) 



2£ 



(z 2 -6 2 (1~ a 2 )) 
(z + l)(l-S 2 ) 



(15) 



For a given initial condition, I and 5 are known, then a is the root of ( |l5| ) that changes with z from its initial value 
continuously. We find that a slightly decreases as the holes get close to each other for most of the initial conditions. 
For example, considering identical holes (S — 0) and setting £ = we find a 2 — (2z + 1) / (z + l) 2 . In Fig 1 whe 
show the change of a as a function of the separation z for 77 = 1/3 with initial value a = 1/3 at z = 40. 
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FIG. 1. Change of a as a function of the separation d for r\ — 1/3 with initial value a — 1/3 at d = 40. 

We shall now consider the torque of the semi-infinite strings on the holes. The strings total angular momenta are 
not finite. Nevertheless, the rate of their change may be computed as follows 



S 2 z 2 
Si zi 

where Sk , k = 1, 2 are the strings spins given by || 

2mia2 d 



m 2 a 2 
mi ai 



(16) 
(17) 



51 — Si- 2 

5 2 = Si- 2 

Sl-2 : 



A 2 
2m 2 cii d 



A 2 

2mim 2 a 2 
A 3 



(m\ — (mi + d) 2 ^j + a\ (mi — (m 2 + d) 2 ^j - (ai — a 2 ) (a\ - a 2 ) 



where A 2 = d 2 — (mi + m 2 ) 2 + (<zi - a 2 ) 2 . That is, we assume that each string exchanges angular momentum with 
its attached black hole. We use a fixed center of mass reference system in which z 2 — zi = d and m±Zi + m 2 z 2 = 
and a further requirement that the holes angular momentum £h = micti + m 2 a 2 is kept constant along the collision 
path. Then the equations (EftO) are equivalent to the algebraic equations 



miS 2 — m 2 Si = 
miai + m 2 a 2 — £h 



(18) 
(19) 



Thus, ai (d) — ai (^initial) + A(d) — A (di„itiai) where A (d) is the solution of (JT8- 19 ) , that is, it is the real root of 
P(A) = 4 (S 2 - 1) A 3 - AI (S + 2S 2 - 3d) A 2 

+ Ul 2 (2(5 + S 2 - 3d) + d (d 2 - 1) (25 - 1) + (d + l) 2 {5 4 - S 2 (d 2 - l) - 2d 3 )] A 

-I ((d +1)(6- l) 2 (S 2 + dS-d 2 + d)+ AI 2 (6 - d)) = 

Here again we set mi + m 2 to one and 6 = m 2 — mi- We find that the holes spins decreases as the holes get close to 
each other. In Fig. 2 we show the change of spins during the collision process. 
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FIG. 2. Change of spins during the collision process. The settings are the same of Figure 1 except for the initial values 
ai = —1/3 and a 2 — 1/6 

The change in the spins are small during most of the collision process, so we neglet the change in the gravitational 
radiation output presented in the Table I. 

In both cases, strut and string's, we see the transfer of angular momentum from the defects to the holes. One then 
expects some exchange of angular momentum between the holes along the collision path. 

VI. CONCLUSIONS 

We developed a semi-analytical approach of treating the head-on collision problem of two Kerr black holes. This 
is based on an exact solution of Einstein's field equations that can be interpreted as the configuration of two Kerr 
black holes with some singularities on the axis from which we obtain the attraction force between the holes. We use 
the force and Newtonian equations of motion to find the dynamics of the holes. We then calculate the amount of 
gravitational radiation released in the collision using the quadrupolc radiation formula. We also discuss the possible 
exchange of the holes spins. 

It is well known that the efficiency of conversion of rest mass into gravitational radiation in the head-on collision 
problem is extremely small for Schwarzschild black holes. We expect higher efficiency for Kerr black holes. If we 
consider two identical spinning holes far apart and compare their initial energies with the final energy of a single Kerr 
hole, assuming hole's angular momentum conservation, the maximum radiation efficiency is 29% for parallel spins 
whereas the maximum ranges between 29% and 50% for antiparallel spins Q . Our results show that the process is 
more efficient for the case of Kerr black holes than for the Schwarzschild black holes, but still small. The Kerr holes 
emit at most 3% of the rest energy fj, 2 /M where /j, and M are the reduced and total mass respectively whereas the 
Schwarzschild holes emit no more than 1% B. Recall that fi/M < 1/4, so that the maximum radiation efficiency we 
found is less than 0.2% of the total rest energy M. 

According to the area theorems for coalescing black holes the antiparallel case can radiate more energy than the 
parallel case |Ti[ ]. Our approach agree with the area theorems. See the first (last) columns of Tabic I for the 
gravitational wave output of parallel (antiparallel) holes in collision. 

We also found that the spin-spin interaction gives a repulsive contribution to the force for both parallel and 
antiparallel spins. In contradistinction, according to Maxwell theory, two magnetic dipoles repel (attract) each other 
and, according to some authors, a spinning test particle is attracted to (repelled from) a large rotating body if they 
are aligned in parallel (antiparallel) Q. 

We then addressed the possible exchange of spins between the holes. With simple and reasonable assumptions we 
conclude that the holes may spin down or spin up depending on the initial conditions and on the constant parameters 
of the configuration. The effect is small for most of the collision path. When the holes are close enough the spins 



may then change up to 100%. It is known that the spin-spin force is of post-post- Newtonian magnitude |lq] . Further 
investigations must be done towards the clarification of this spin exchange effect. 
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e J -> 


1 


3/4 


1/2 


1/4 





i/A 


i/2 


3i/4 







9.0948 


13.1064 


16.3233 


18.3563 


19.0476 


19.7450 


21.8679 


25.4811 


30.6095 


-0.1 


7.4681 


10.9436 


13.7352 


15.5004 


16.1008 


16.7065 


18.5506 


21.6896 


26.1450 


-0.2 


5.4081 


8.2322 


10.5088 


11.9502 


12.4408 


12.9357 


14.4430 


17.0095 


20.6521 


-0.3 


3.0990 


5.0773 


6.6864 


7.7086 


8.0569 


8.4084 


9.4797 


11.3052 


13.8951 


-0.4 


1.0577 


1.9689 


2.7258 


3.2101 


3.3755 


3.5427 


4.0526 


4.9225 


6.1549 



TABLE I. In the first column are the values of e. In the first line are the values of 0. The other entries of the table give the 
total amount of gravitational radiation emmitted during the collision process. 
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